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Abstract 

In the framework of nonrelativistic QCD, we compute the leading double- logarith- 
mic corrections of order \T?{l/as) to the heavy-quark-antiquark bound-state wave 
function at the origin, which determines the production and annihilation rates of 
heavy quarkonia. The phenomenological implications for the top-antitop and T 
systems are discussed. 
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1 Introduction 



Nonrelativistic quantum chromodynamics (NRQCD) is a powerful tool for the in- 

vestigation of heavy-quark threshold dynamics. Recent developments of the NRQCD 
effective-theory approach [^^,|6|,|7|,|8|j9|, p!0| . [TT| , p!2| , p!3| , p!^ led to the complete next-to-next-to- 
leading-orderQ (NNLO) description of the production of heavy quark-antiquark pairs at 
threshold. Important applications include T sum rules and toponium phenomenology 
P,|T|,|l3,|ll,|l9|,||,|l],|2|,|||l|5l|6|,|3. A review of the recent progress in the perturba- 



tive study of heavy quark-antiquark systems may be found, for example, in Ref. In 
view of the surprising significance of the NNLO corrections, it appears indispensable to 
also gain control over the next-to-next-to-next-to-leading order (N^LO) in order to im- 
prove the reliability of the theoretical predictions and our understanding of the structure 
and the peculiarities of the threshold expansion. 

Some specific classes of N^LO corrections were analyzed in literature. The one-loop 
renormalization of the 1/^^ operators was obtained in Refs. [^,^. The retardation 
effects arising from the emission and absorption of virtual ultrasoft gluons by the heavy 
quarks were studied in Ref. 0]. The leading logarithmic O {a^\n{l/as)) corrections to 
the heavy-quark bound-state energies En, where n is the principal quantum number, were 
obtained in Ref. . 



In this paper, we take the next step in this direction and investigate a particular 
class of N^LO corrections, namely the leading logarithmic O (^a^ \n'^{l/as)^ corrections to 
the wave functions at the origin ipn{0) of the heavy quark-antiquark bound states which 
are not generated by the renormalization group (RG). As is well known, ipn{0) are key 
parameters in the analysis of the creation and annihilation of heavy quarkonia. The origin 
of these logarithmic corrections is the presence of several scales in the threshold problem. 
In fact, ln{l/ as) appears as a logarithm of a ratio of scales. These corrections are related 
to the anomalous dimensions of the operators in the effective Hamiltonian. They can 
be found by analyzing the divergences of the effective theory or by direct inspection of 
the regions of the logarithmic integration. We shall verify that both methods lead to the 
same result. As a by-product of our analysis, we shall also reproduce the O {a^\n{l/as)) 
corrections to En recently obtained in Ref. 

In N^LO, there are in general also logarithmic corrections of the form ln"^{fi/asmq) 
{m = 1,2,3), where agmg represents the soft or potential scales (see Section 2). These 
corrections may be directly extracted from the NNLO result via the RG equation, and 
they may be resummed by an appropriate choice of the normalization point, /i ~ asm^. 
We note in passing that, in the MS scheme, the optimal choice in NNLO is /i ~ few units x 



This paper is organized as follows. In Section 2, we recall the potential-NRQCD 
(pNRQCD) formalism, from which En and ipn{0) can be extracted. In Section 3, we 
evaluate the non-RG O (^a^ ln^(l/as)) and O (a^ln(l/as)) corrections to ipni^) and En, 

^In the NRQCD effective theory, there are two expansion parameters, the strong couphng constant 
as and the heavy-quark velocity /3, and the perturbative order of some correction is determined by their 
total power. E.g., terms of 0{al), 0{asP), and 0(/3^) contribute at NNLO. 
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respectively, using the effective-theory approach. In Section 4, we repeat this calculation 
in the conventional approach, by inspecting the logarithmically divergent phase-space 
integrals. In Section 5, we present a numerical analysis and discuss phenomeno logical 
implications of our results. Section 6 contains our conclusions. 



2 Vacuum-polarization function near threshold in 
NRQCD 

We investigate the near-threshold behavior of the vacuum-polarization function n(g^) of 
a heavy-quark vector current = q'j^q, 

{q,q. - 9,uq') Il{q') = i J rf^xe^''-^(0|Tj^(a;)j.(0)|0). (1) 

Its imaginary part is related to the normalized cross section of qq production in e+e~ 
annihilation at energy s = q^, 

= ;f^'-f) (2) 
cr(e+e ) 

by 

R{s) = 127rQjlmn(s + ze), (3) 

where Qq is the fractional charge of quark q. 

Near threshold, t he heavy qu arks are nonrelativistic, so that one may consider the 
quark velocity f3 = ^1 — Am'^/s as a small parameter. An expansion in jS may be per- 
formed directly in the Lagrangian of QCD by using the framework of effective field theory. 



In the threshold problem, there are four different scales [12]: (i) the hard scale (energy 
and momentum scale like mq); (ii) the soft scale (energy and momentum scale like (3mq); 

(iii) the potential scale (energy scales like P'^niq, while momentum scales like Prriq)] and 

(iv) the ultrasoft scale (energy and momentum scale like jS'^rriq). The ultrasoft scale is 
only relevant for gluons. By integrating out the hard scale of QCD, one arrives at the ef- 
fective theory of NRQCD l^l^jQ- If one also integrates out the soft scale and the potential 
gluons, one obtains the effective theory of pNRQCD, which contains potential quarks and 



ultrasoft gluons as active particles |T^. The dynamics of the quarks is governed by the 



effective, nonrelativistic Schrodinger equation and by their interaction with the ultrasoft 
gluons. To get a regular perturbative expansion within pNRQCD, this interaction should 
be expanded in multipoles. The corrections from harder scales are contained in the Wilson 
coefficients, leading to an expansion in a^, as well as in the higher- dimensional operators 
of the nonrelativistic Hamiltonian, corresponding to an expansion in l/rUq or f3. 

The nonrelativistic expansion in as and (3 provides us with the following representation 
of the heavy-quark vacuum-polarization function near threshold: 

N 

n(i?) = ^C,(a,)G(0,0,E) + ..., (4) 
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where E — ^/s — 2m q is the qq energy counted from the threshold, Ch{(y.s) is the square 
of the hard matching coefficient of the nonrclativistic vector current, and the elhpsis 
stands for the higher-order terms in j3. G'(x, y, E) is the nonrclativistic Green function, 
which sums up the (ag/f^)^ terms singular near the threshold. It is determined by the 
Schrodinger equation which describes the the propagation of the nonrclativistic quark- 
antiquark pair in pNRQCD, 

(7i-E)G(x,y,^) = 5(3)(x-y), (5) 

where Ti. is the nonrclativistic Hamiltonian defined by 

^ + V{x) + ..., V{x) = Vc{x) + ... . (6) 

mg 

Here, Vc{x) = —Cpag/x is the Coulomb potential, Cp = 4/3 is the eigenvalue of the 
quadratic Casimir operator of the fundamental representation of the colour group, x = |x|, 
and the ellipses stand for the higher-order terms in ag and l/mg. The Green function has 
the spectral representation 

G{^,y,E)-}^^ i^n-i^ +io kymg - E ' 

where iprn a-nd ip^ are the wave functions of the qq bound and continuum states, respec- 
tively. 

Below the threshold, the (perturbative) vacuum-polarization function of a stable heavy 
quark is essentially determined by the bound-state parameters. For the leading-order 
Goulomb (G) Green function, the energy levels and wave functions at the origin read 



c _ '^s 



E": = 



m„n^ 



2 



A! 



where As = asCF'mq/2. We are interested in the corrections to ■^^^^(O) . Note that, for 
the study of bound-state parameters, we have /3 ~ cts, so that we are only dealing with 
one expansion parameter. 



3 Non-RG leading logarithmic corrections from the 
effective- theory approach 

The NNLO non-RG leading logarithmic corrections to the wave functions at the origin 
are generated by the following operators in the effective nonrclativistic Hamiltonian: 

A'„ = + ^ ij + (l + Is') ^%^i(x) - -il. (9) 

2mqX^ 2ml ^ ^ 3 / ^ ' Am^ ^ ' 
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where Ca = 3 is the eigenvalue of the quadratic Casimir operator of the adjoint repre- 
sentation of the colour group, S represents the spin of the quark-antiquark system, and 
{., .} denotes the anticommutator. The first term of Eq. {^) is the non-Abelian potential 
||3^ , |33| , |3^ , and the rest is the standard Breit potential and the kinetic-energy correction. 
Note that the representation of the nonrelativistic Hamiltonian is not unique and 
can be related to other representations found in literature by the use of the equations 
of motion. However, the corresponding corrections to the Green function [|T8|JT9|j2^ are 
independent of the specific representation. In the vicinity of the nth bound-state pole, 
the resulting correction to the Green function reads [^] 



AG(0,0,E)|^^„c 



EC-E ml 



A-^S{S + l) 



Cf + 2Ca\Gc{^,Q,E) + 



(10) 

where Gcl^, y,£') is the Coulomb Green function. The spin-independent term propor- 
tional to Cf in Eq. (|10D comes about as (4 — 1-1- 1)Cfi where the contributions arise from 
the anticommutator, 5-function, and l/m^ terms of Eq. (^. Writing 



|V'n(0)P=|^^(0)f (l + A^^ 



we have 



AV'^(O) 



A--S{S + l] 



Cf + 2Ca\Gc{^,Q,E^) + ... . 



fill 



(12) 



These corrections are singular, since the Coulomb Green function at the origin is ultravio- 
let (UV) divergent. In dimensional regularization, with (i = 4 — 2e space-time dimensions, 
it takes the form 



Gc;(0,0,A;) 



1 

Ye 



k 



2 



+ 



rUgX 



27r V2e k 



+ 



(13) 



where k"^ 



-TUqE and the ellipsis stands for the nonlogarithmic contribution. The 



pole term in Eq. (|TB|) is canceled by the O {a1) infrared (IR) pole of the hard matching 
coefficient Ch{ois) [0,11^. Thus, the scale /x in the logarithm is to be identified with m, 



and the sought correction reads P^pT | 



AV^(O) = CfciI 



2--S{S + l) 



Cf + Ca\ In 



a. 



(14) 



The residual operators in the NNLO effective nonrelativistic Hamiltonian which are not 
contained in Eq. (|^) correspond to the purely perturbative corrections to the static 
Coulomb potential. The corresponding corrections to the wave functions at the origin 
T9| , |2l| , p6| contain RG logarithms of the form ln™(/i/asmg) (m = 1,2), which vanish 
This also holds in N^LO for the RG logarithms of the form asTnq) 



a.m. 



for /i 

(m = 1, 2, 3) because the ultrasoft effects enter the stage only in N'^LO, so that the cor- 
responding running of the strong coupling constant at the ultrasoft scale only becomes 
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relevant in N^LO. An important point here is that, starting from NNLO, the hard match- 
ing coefficient Chicks) receives a nonvanishing anomalous dimension. Therefore, starting 
from N^LO, not only the running of should be taken into account, but also the effective- 
theory RG should be used for the evolution of the hard matching coefficient from = ruq 
down to /i = agmq [p3| . 

At N^LO, the non-RG leading logarithmic corrections to the wave functions at the 
origin are produced by the one-loop renormalization of the operators in Eq. (P). In 
dimensional regularization, the pole part of the correction is 




4 2 



+ 



2„ /17 



-Cf + 

3 V 3 3 



Ca 



nigX 



2 + 3 



X 



16 8 \ TTQs 



5(x) 



(15) 



where the ffist three terms contained within the parentheses represent the IR divergence, 
while the fourth one embodies the UV divergence of the potential. The IR poles are 
canceled by the ultrasoft contribution, with the characteristic scale almq, and may be 



read off from Refs. ||30| , ^ , while the UV poles are canceled by the IR poles of the hard 
coefficients and may be extracted from Refs. 0,^. Evaluating the corrections to the 
wave functions at the origin due to Eq. ([l^) in the same way as Eq. (0) was obtained, 
we find the N^LO non-RG leading logarithmic corrections to the wave functions at the 
origin to be 



A>^(0) 



TT 



41 7 



CfCa + -Cl In' 



(16) 



In the derivation of this result, the factor {fi/kY^ in Eq. (0) needed to be expanded up 
to O(e^). The contributions to Eq. (0) related to the IR and UV poles of Eq. ( |T5| ) are 
of opposite signs. This may be understood by observing that the IR poles introduce a 
factor ln(£'^/As) ~ Ina^, while the UV poles contribute a factor ln(mq/As) ~ ln(l/a;s). 
The Abelian Cp term in Eq. ([I6|) agrees with the QED result of Ref . |^ . In contrast to 
the QED case, the leading logarithmic QCD corrections are spin dependent. 

In the remainder of this section, we revisit the leading logarithmic corrections to the 
bound-state energy levels, which may be obtained in a way similar to the case of the wave 
functions at the origin. These corrections start from N^LO, since there are no relevant 
singularities in NNLO. In addition to Eq. (|15D, we now have to take into account the 
IR-singular contribution to the Coulomb potential [§^j3l 

1 CpClai 
2e 12x ' 



which is canceled by the ultrasoft contribution 
butions to the wave functions. Writing 



(17) 

Obviously, this term gives no contri- 



(1^ 
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we find 



TT in 



41 7 2 
.on on 3n"^ 



2 4 2-^3 

CpCA + i^—CfC^ + -C^ 
dn 



1 

— ( 
6 



X In — , 

as 



(19) 



in agreement with the result for / = of Ref. [31 



4 Non-RG leading logarithmic corrections from the 
traditional approach 

An ahernative method of finding the leading logarithmic corrections is to directly in- 
spect the regions of logarithmic integration [^,^,^,^. Let us first consider the NNLO 
corrections to the wave functions. Substituting the continuum part of the spectral repre- 
sentation (1^) into Eq. (|12D, we obtain 



AV^(O) 



m„ 



A--S{S + 1] 



(27r)3 k'^/mq - E„ 



(20) 



The integral over k logarithmically diverges at large momentum and should be cut at scale 
rriq, where the nonrelativistic approximation becomes inapplicable. At low momentum, 
the logarithmic integration over k is effectively cut at the Coulomb scale agmg. Inserting 
in Eq. (|20|) the well-known expression for the Coulomb wave function at the origin, 



27rA, 



1 



1 + 



TT 



A. 



+ ... 



(21) 



k 1 — exp (— 27rAs//c) ' k 
where only the second term is relevant for our purposes, we find, with logarithmic accuracy. 



A'<(0) = C^ai 



2 --5(5 + 1) 



™9 dk 



Ctsiriq 



k 



+ 



(22) 



Performing in Eq. (^2]) the integration over k, we recover Eq. (0). 

The N^LO leading (double) logarithmic corrections to the wave function at the origin 
are of the form 



A>n(0) 



'Cpal < 

-'.CI 



°-CI + ACaCf + ^CI 
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S{S+1)]CfCa 



1 dk r'^ dk' 

asTUq k J — E k' 

•^1 dk dk' 

?.mo k Jk k' 



(23) 



Here k is the momentum of the potential heavy quark, which can be as small as agrrig, 
while k' is the momentum of the virtual gluon, which can even be ultrasoft, of order 
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arm, 



The integrals over k have the same origin as those in Eq. (|2^). The integrals over 
k' represent the logarithmic corrections to the potential due to the virtual-gluon exchanges 
Il5| j29| . |30| ■ In the first integral on the right-hand side of Eq. (p3D, the quark momentum k 



plays the role of an UV cutoff, and the quark energy —E = k /niq acts as an IR cutoff, 



so that this contribution corresponds to the IR poles of Eq. (pTq). In the second integral, 
the momentum k acts as an IR cutoff, and this contribution corresponds to the UV poles 
of Eq. (|T5|) . Integrating Eq. ( p3|) over k' and k, we arrive at Eq. ([161) . A similar analysis 
for QED bound states may be found in Ref. [pS] . 

The N^LO leading (single) logarithmic corrections to the energy levels may be obtained 
by the method introduced in Refs. p8| , p9t^ , where the regions of the virtual-photon 



momentum which lead to logarithmic contributions have been studied in the QED case. 
In this way, we obtain 




+ 



-cl 



17 



31 
Qn 
7 



2 



'^^F + ( ^ ) ^F^^ + ^CfC\ 

in \bn 671'^ J 6n 



Qn 



5(5 + i; 




dk' 



(24) 



where again the first (second) integral corresponds to the IR (UV) poles of Eq. ([T3|). After 
integration, we recover Eq. (plOl). 



5 Phenomenological applications 

Our results affect two processes of primary interest, namely the threshold production of 
top and bottom quark-antiquark pairs. The relatively large width Ff of the top quark 
serves as an efficient IR cutoff for long-distance effects. Because the relevant scale ^/rnt^l 
is much larger than the asymptotic scale parameter Aqcd? but comparable to the Coulomb 
scale, the cross section in the threshold region may be described by the NRQCD pertur- 
bation expansion if singular Coulomb effects are properly taken into account ET|. In 



order to analyze the significance of the N^LO leading logarithmic corrections to the cross 



section, we start from the NNLO calculation of Ref. [261 and add the contributions from 



Eqs. (|T^) and (|1^). In Fig. 1, the normalized cross section R thus obtained is compared 
with the pure NNLO result. The input parameters are taken to be as{Mz) = 0.118, 

= 175 GeV, and F^ = 1.43 GeV. The soft renormalization scale is determined from the 
condition fig = 2as{^is)^ti and hard renormalization scale is chosen to be /i^ = m^. The 
effect of the N'^LO leading logarithms is twofold. The normalization of the cross section 
is reduced by about 7% around the IS peak and below, and the energy gap between the 
IS peak and the nominal threshold is decreased by roughly 10%. 

In the case of bottom quark-antiquark production, the nonperturbative effects are 
much more significant, and one is led to use the sum rule approach |^^,^ to get them 
under control. Specifically, appealing quark-hadron duality, one matches the theoretical 
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results for the moments of the spectral density, 



Mn 



n 



(4m 



s=0 



2\n 



ds- 



.71+1 



(25) 



with their experimental counterparts evaluated from the expressions in the second line of 
Eq. (^). For large n, the moments are saturated by the near-threshold region. Then, 
the main contribution to the experimental moments comes from the T resonances, which 
are measured with high precision. On the other hand, for n of 0(l/a^), the Coulomb 
effects should be properly taken into account on the theoretical side. In order to analyze 
the N^LO leading logarithmic corrections to the T sum rules, we upgrade the NNLO 
result of Ref. ^T] by including Eqs. ([I6| ) and (pJ)). We fix the strong coupling constant by 
as{Mz) = 0.118 and focus on the determination of the bottom-quark mass. At present, 
this appears to be the most interesting application of Eq. (pSl) . We find that the inclusion 
of the N^LO leading logarithms in the sum rules leads to a reduction of the extracted 
mass value by approximately 30 MeV for moderate values of n, 5 < n < 15, and if 
the soft renormalization point Hs is chosen to be fig = 4as(/is)mfe. This result does not 
depend on whether the energy denominators of the Green function are expanded around 



the Coulomb values or not (see Ref. for details) and on which mass parameter, pole or 
MS mass, is considered. On the other hand, the result essentially depends on /i^ because 
the /is dependence of is compensated only by higher-order terms. For example, for 
fig = 2as(/is)mb, the correction to the mass parameter reaches —70 MeV. However, the 



perturbative result can hardly be trusted at such a low renormalization point [21, 2(: 



For completeness, we also present the numerical corrections to the parameters of the 
IS T resonance. For = 4.8 GeV and fig = 4as(/is)mfe, the wave-function correction is 
— 19%, and the correction to the binding energy is 25 MeV. 



6 Discussion and conclusions 

We studied a special class of NNLO and N^LO corrections to the key parameters of heavy 
quark-antiquark bound states, namely those which are enhanced by a maximum power of 
ln(l/a;s) and are not generated by the RG. By contrast, the RG logarithms are well known 
and may be resummed by an appropriate scale choice. Such non-RG leading logarithmic 
corrections first arise for the wave functions at the origin in NNLO |T9| , ^ and for the 
energy levels in N^LO [^. Specifically, they are of the forms ln(l/as) and ln(l/Q;s), 
respectively. We confirmed these results and completed the knowledge of non-RG leading 
logarithmic corrections in N^LO by computing the O (^af ln^(l/as)^ corrections to the 
wave functions at the origin. 

We applied our result to top and bottom quark-antiquark production at threshold. In 
the case of top, the resulting correction to the production cross section near the IS" peak 
reaches 10%. In the case of bottom, the mass extracted from the T sum rules is shifted by 
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approximately —30 MeV. The latter value is comparable to the uncertainty of ±60 MeV, 
which is usually assigned to this kind of bottom-quark mass determination on the basis 
of the renormalization scale dependence of the strong coupling constant in the NNLO 
result. In fact, the non-RG leading logarithmic contributions are about two times smaller 
than the contribution from the N'^LO RG logarithms, which may be estimated from the 



renormalization scale dependence of the NNLO result ||21| , |26| . Thus, the resummation of 



the RG logarithms |^ seems to be very justified. At the same time, the scale of the 
corrections is close to the estimate given in Ref. on the basis of the analysis of the 
ultrasoft contributions. 

Although, in contrast to QED, ln(l/Q;s) is not a big number, especially for the case of 
bottom, the leading logarithmic terms can be considered as typical representatives of the 
N^LO corrections. For comparison, at NNLO, the leading logarithmic term accounts for 
approximately one half (third) of the total correction to the n = wave function at the 
origin in the case of top (bottom). Obviously, the N^LO corrections are comparable to 
the NNLO ones and reach 10% in magnitude, even in the case of top, where ~ 1/10. 
This tells us that the NRQCD threshold expansion is not a fast convergent series for the 
physical value of the strong coupling constant. 

A final comment refers to the resummation of the non-RG leading logarithmic cor- 
rections to the wave functions at the origin. A part of the non-RG leading logarithmic 
corrections to the heavy-quark threshold cross section was resummed in Ref. |^3| by using 
the RG equation of the effective theory for the evolution of the hard matching coefficient 
Chi^s) from scale nig down to scale Prrig. This evolution equation is obtained by studying 
the dependence of Chids) on the scale fi, which cancels the /i dependence of Eqs. (|T2|) and 
(p!3|). This effectively sums up the higher-order NRQCD corrections due to the tree- level 
operators of Eq. (^). The terms resummed in this way are of the form as{ois ln(l//3))^"'~^, 
with n = 1, 2, . . ., i.e. they include only even powers of a^. For the bound-state parame- 
ters, we have /3 ~ a^, and Eq. ( |T6| ) gives the first term of this series. In order to resum all 
correction of the form as(as ln(l/as))"', it is necessary to add the terms with odd powers 
of as, of the form ^^(as ln(l/as))^"', which are generated by Eq. (plSj) . For this end, one 
has to take into account not only the evolution of Ch{c(s), but also the evolution of the 
potential (|^) from the hard scale down to the ultrasoft scale. Numerically, however, the 
effect of the resummation may not be essential for phenomenological applications. 
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Fig. 1. Normalized cross section R{E) of e^e~ — > tt in NNLO (dotted line) and with the 
leading logarithmic N^LO corrections included (solid line), as a function of the centre-of- 
mass energy E counted from the nominal threshold at 2mt. 
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